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Abstract — Precise tracking of a randomly varying optical 
phase is key to metrology, with applications in optical commu- 
nication. Continuous phase estimation is known to be superior 
in accuracy as compared to static estimation. The underlying 
parameters in the signal model or the measurement itself 
are, however, prone to changes owing to unavoidable external 
noises or apparatus imperfections. The estimation process is, 
thus, desired to be made robust to uncertainties in these 
parameters. Here, homodyne phase estimations of coherent 
and squeezed states of light, evolving continuously under the 
influence of a resonant noise process, are made robust to 
parameter uncertainties using a robust fixed-interval smoother, 
designed for uncertain systems satisfying a certain integral 
quadratic constraint. 

I. INTRODUCTION 

The quantum phase estimation [1] is central to various 
fields such as metrology [2], quantum computation [3], 
communication [4], [5] and quantum cryptography [6], A 
fundamental bound on the estimation accuracy is imposed 
by Heisenberg's uncertainty principle [7], that, for example, 
limits gravitational wave detection [8]. A less fundamental 
yet practically relevant limit, called the standard quantum 
limit (SQL), is the minimum level of quantum noise obtained 
using standard approaches without real-time feedback. The 
SQL sets an important benchmark for the quality of an 
estimation process and calls for devising quantum strategies 
that can beat it. 

Real-time feedback in homodyne estimation of a static 
unknown phase can yield mean-square errors lower than the 
SQL [9], [10], [11], [12] or even reaching the theoretical 
limit [13]. However, it is experimentally more relevant to 
precisely estimate a phase, that is continuously varying under 
the influence of an unmeasured classical noise process [14], 
[15], [16], [17]. A classical process coupled dynamically to 
a quantum system under continuous measurement may be 
estimated in various ways: prediction, filtering or smoothing 
[18]. In particular, smoothing uses both past and future 
measurements to yield a more accurate estimate than filtering 
alone, that uses only past measurements. Since smoothing is 
non-causal, it cannot be used in real-time, but is used for 
offline data processing or with a delay with respect to the 
estimation time. 

The fixed-interval smoothing problem [19], [20], [21] 
considers measurements over a fixed time-interval T, where 
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the estimation time is t < T. One solution, the Mayne-Fraser 
two-filter smoother [22], [23], [24], consists of a forward- 
time Kalman filter and also a backward-time Kalman filter 
called an "information filter" [25] and it combines the two 
estimates to yield the optimal smoothed estimate. The infor- 
mation filter and the smoother were combined into a single 
backward smoother by Rauch, Tung and Striebel (RTS) [26], 
[27]. 

In Refs. [16], [17] the signal phase to be estimated is 
allowed to evolve under the influence of an unmeasured 
continuous-time Ornstein-Uhlenbeck (OU) noise process. 
However, the estimation process heavily relies on the un- 
derlying parameters being precisely known, which is practi- 
cally not feasible due to unavoidable external noises and/or 
apparatus imperfections. It is, therefore, desired to make the 
estimation process robust to uncertainties in these parameters 
to improve the quality of the phase measurement. In Ref. 
[28] the authors have shown that a robust guaranteed cost 
filter [29] yields better accuracy in the phase estimate than 
an optimal Kalman filter with parameter uncertainty in the 
phase being measured. The authors have also demonstrated 
improvement in continuous phase estimation with robust 
fixed-interval smoothing [30] as compared to optimal RTS 
smoothing for a coherent state [31] and a squeezed state [32] 
of light. 

These works, however, considered a simplistic Ornstein- 
Uhlenbeck noise process modulating the signal phase to be 
estimated, and the kind of noises that in practice corrupt 
the signal are more complicated than an OU process. The 
authors have illustrated in Ref. [33] a guaranteed cost robust 
filter designed for a more complicated and practically rele- 
vant second-order resonant noise process. Here, we design 
a robust fixed-interval smoother [30], meant for uncertain 
systems admitting a certain integral quadratic constraint, for 
coherent light beam phase-modulated by such a resonant 
noise process. We compare the performance of the robust 
and RTS smoothers for the uncertain system with respect to 
the standard quantum limit. We also extend this by designing 
such a robust fixed-interval smoother for a phase-squeezed 
light beam and compare its behaviour with a corresponding 
RTS smoother in relation to the coherent state limit (CSL). 
The coherent state limit determines the minimum estimation 
error that can be obtained for a coherent state of light. 

II. RESONANT NOISE PROCESS 

The resonant noise process under consideration in this 
paper is typically generated by a piezo-electric transducer 
(PZT) driven by an input white noise. The simplified transfer 



function of a typical PZT is given by: 
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where k is the gain, £ is the damping factor, uj r is the 
resonant frequency (rad/s) and u is a zero-mean white 
Gaussian noise with unity amplitude. 

We use the following values for the parameters above: 
k = 9 x 10 4 , C = 0.01 and uj r = 6.283 x 10 3 rad/s (1 kHz). 
Fig. 03 shows the Bode plot of the transfer function ([T}. 

Let X\ := cf> and x-i :— (f>. A state-space realization of the 
transfer function (Q3 is: 



x = Ax + Gv, 



(2) 
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COHERENT STATE: OPTIMAL SMOOTHER 

The homodyne "smoothed" phase estimation of a weak 
coherent state of light is optimal for given values of the 
parameters when using an offline RTS estimator alongwith 
a Kalman filter in the feedback loop to adjust the phase of 
the local oscillator. Under a linearization approximation, the 
homodyne photocurrent from the phase estimation system is 
given by [16]: 



I(t)dt = 2\a\[<p(t) - 4>(t)]dt + dW(t), 



(3) 



where \a\ is the amplitude of the coherent state with photon 
flux given by J\f := \a\ 2 , <j> is the intermediate phase 
estimate, and W(t) is a Wiener process arising from the 
quantum vacuum fluctuations. 

The measurement model is, thus, [33]: 



Hx + Jw, 



(4) 



where w := is a zero-mean Gaussian white noise with 
unity amplitude and H := I 2\a\ J and J := 1. 

Rewriting the equations for the system under considera- 
tion, we get: 

Process model: x = Ax + Gv, 
Measurement model: 8 = Hx + Jw. 



(5) 
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Fig. 1. Bode plot of the resonant noise process transfer function. 



E[v(t)v{r)] = N8{t - 
E[w(t)w(T)} = S6{t- 
E[v(t)w{r)] = 0. 
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Since v and w are unity amplitude white noise processes, 
both N and S are unity (scalars). 

A. Forward Kalman Filter 

The continuous-time algebraic Riccati equation to be 
solved to construct the steady-state forward Kalman filter 
is [34]: 

AP f + P f A T + GNG T -PfH T (JSJ T ) ^HPf = 0. (6) 
The forward Kalman filter equation is [34]: 

£ = (A - K f H)x + K f Hx + K f Jw, (7) 
where K f := P f H T (JSJ T ) _1 is the Kalman gain. 

B. Backward Kalman Filter 

The continuous-time algebraic Riccati equation to be 
solved to construct the steady-state backward Kalman filter 
is [27]: 



AP b - P b A T + GNG T P b H T (JSJ T ) _1 HP b 

The backward Kalman filter equation is [27]: 

£ = ( A K b H)x + K b Hx + K b Jw, 

where K b := P b H T (JSJ T ) _1 is the Kalman gain. 

C. Smoother Error 

The smoother error covariance matrix is obtained as 



0. (8) 



(9) 



(10) 



since the forward and backward estimates of the optimal 
(Kalman) smoother are independent [27]. 

Using k = 9 x 10 4 , C = 0.01 and uj r = 6.283 x 10 3 rad/s 
(1 kHz) as in section Ull and |a| = 5 X 10 2 , we get: 



3.7998657 x W J 








(11) 



3.7344282 x 10 5 
IV. COHERENT STATE: ROBUST SMOOTHER 

In this section, we make our smoother robust to uncertain- 
ties in the parameters underlying the system matrix A using 
the robust fixed-interval smoothing approach from Ref. [30] . 

We introduce uncertainty in A as follows: 



A^ A 
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where A := [ 8% 82 1 is an uncertain parameter satisfying 
||A|| < 1 that implies 8\ + 8% < 1. < fix < 1, < /i 2 < 1 
are parameters determining the levels of uncertainty. 

Thus, the process and measurement models of (0 take the 
form: 



Process model: x = (A + GAK)x + Gv, 
Measurement model: 6 = Hx + Jw, 



(12) 



where K := 











The Integral Quadratic Constraint of Eq. (2.4) in Ref. [30] 
for our case is: 



(uf+v A )dt < 1 + 

JO 



l dt, 



(13) 



where z = Kx is the uncertainty output, and w — AKx + v 
and v = w are the uncertainty inputs. Thus, we would have 
Q = 1 and R = 1 in our case. 

The steady-state forward Riccati equation, as obtained 
from Eq. (5.1) in Ref. [30] for our case, is: 

YA + A T Y + YGQ 1 G T Y + K T K H T RH = 0. (14) 

The steady-state backward Riccati equation, as obtained 
from Eq. (5.2) in Ref. [30] for our case, is: 

ZA + A T Z ZGQ 1 G T Z K T K + H T RH = 0. (15) 

Next, Eq. (5.3) in Ref. [30] yields: 

fj = -(A + GQ _1 G T Y) T ;; + H T R6>. (16) 

Likewise, Eq. (5.4) in Ref. [30] for reverse-time yields: 

£ = (A - GQ _1 G T Z) T £ + H T R0. (17) 

The forward filter is simply the centre of the ellipse of 
Eq. (3.3) in Ref. [30]: Xf = Y _1 r;. Likewise, the backward 
filter is: x b = Tr 1 ^. 

The robust smoother for the uncertain system would, then, 
be the centre of the ellipse of Eq. (5.5) in Ref. [30]: 



x=(Y + Z)- 1 ( 7? -e). 



(18) 



V. COHERENT STATE: COMPARISON FOR THE 
ROBUST AND KALMAN SMOOTHERS 

A. Lyapunov Method 

1) Forward Filter: We augment the system given by 
( fT2l with the forward Kalman filter (0 and represent the 
augmented system by the state-space model: 



= A: 



Bw, 



(19) 



where 

x := 

Thus, we have 
A = 



x 

X 



A + GAK 
K f H 



and 
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For the continuous-time state-space model ([T9T >. the 
steady-state state covariance matrix Pf s is obtained by solv- 
ing the Lyapunov equation: 

-r-T =-=rT 



AP fs +P fs A +BB 
where Pf s is the symmetric matrix 



fs 



£(xx' 



S 



0, 



M f 
N f 



(20) 



The state estimation error can be written as: 

ef := x — x = [1 — l]x, 

which is mean zero since all of the quantities determining 
ef are mean zero. 

The error covariance matrix is then given as: 



n f : = £(e f ef J ) = [1 - l]£(xx 



[1-1] 
S - M f 



£ M f 
M f T N f 

- M f T + N f . 



(21) 



2) Backward Filter: The augmented system state-space 
model $1% for the backward Kalman filter (O would have: 



A = 



A + GAK 
K b H 





A K b H 



G 
K b J 



We then solve the Lyapunov equation d20l . with Pf s 
replaced by 



bs 



£(xx n 



£ 
Ib J 



M V T 



M b 
N b 



for the backward filter. 

The error covariance matrix is, thus: 

n b := E(e h e h T ) = £ - M b - M t 



Nv 



(22) 



3) Cross-Correlation Term: The forward and backward 
estimates are not independent and are correlated in this case, 
unlike in section [TlTl The cross-correlation term is [31]: 



nfb := £(e f e b 



M f 



M b + a£/3, (23) 



where a := M/S" 1 and j3 := £- x M b [21]. 

4) Smoother Error: The smoother error covariance is 
[31]: 



n 



n 



fb 



(24) 



n f + n b 2n fb 

Since we are mainly interested in estimating x\ = (f>, the 
estimation error covariance of interest is a 2 = 11(1, 1). 

B. Standard Quantum Limit 

The standard quantum limit is set by the minimum phase 
estimation error that can be obtained using perfect hetero- 
dyne scheme. The heterodyne scheme of measurement incurs 
the same noise penalty as dual-homodyne scheme [16]. A 
schematic for the dual-homodyne scheme is depicted in 
Fig. An electro-optic modulator (EOM), driven by the 
resonant noise source, phase-modulates the input signal. The 
modulated signal then splits through a 50 — 50 beamsplitter 
into two arms each with a homodyne detector (HD1 and 
HD2, respectively, with the local oscillator phase of HD1 
7r/2 out of phase with that of HD2). The arctan of the ratio 
of the output signals of the two arms is then fed into a low- 
pass filter (LPF). The phase estimation error would then be 
minimum if this LPF is an optimal Kalman filter. 



The output signals of the two arms are: 
h = —j= (2\a\ sin <f> + m + rt 2 ) , 

h = —7= (2|a| cos^ + n 3 - n 4 ) , 
V2 

where 712 and 77,4 are the noises arising from the vacuum en- 
tering the empty port of the input beamsplitter corresponding 
to the two arms, respectively, and n\ and 77,3 are measurement 
noises of the two homodyne detectors, respectively. All these 
noises are assumed to be zero-mean white Gaussian. 
The output of the arctan block is: 

2\a\ sin 6 + m + 77,2 



•d = arctan 



(25) 



2\a\ cos <f> + 77,3 — ri4 / 

Assuming the input noises are small, a Taylor series 
expansion upto first-order terms of the right-hand side yields: 

1 1 



2\a 



■ni 



2\a 



■n,2- 



(26) 



Expressing this equation in terms of x, we get the mea- 
surement model as: 



■d — Hx + J 



(27) 



1 1 

2\a\ 2\a\ 

of the optimal steady-state 



where H = [ 1 ] and J = 
The error covariance matrix 
Kalman filter for the process given by (fJJ and the mea 
surement given by ( l2~7| > is obtained by solving an algebraic 
Riccati equation of the form (O for Pf . The error covariance 
of interest (i.e. that in estimating x\ = <fr) is then a 2 — 
Pf(l,l). 

C. Comparison of Estimation Errors 

The estimation mean-square errors may be calculated, as 
described in section fV^Al for the Kalman smoother, and like- 
wise for the robust smoother, as a function of the uncertain 
parameter Si with 82 — 0. Here, we use the nominal values 
of the parameters and choose different values for /ii and set 
[12 = 0. These values were used to generate plots of the 
errors versus <5i to compare the performance of the robust 
smoother and the Kalman smoother for the uncertain system 
in comparison to the SQL. The SQL value is obtained by 
designing, as described in section IV-BI a different Kalman 
filter for each value of the uncertain parameter. Figs. |3] |4] an d 
[5] show these plots for [i\ = 0.25, 0.5 and 0.75, respectively. 

Clearly, the Kalman smoother behaves better than the 
robust smoother when 5i — 0, as expected. However, as <5i 



Resonant noise 



Signal 



EOM 



a I e 





HD1 












HD2 







Sv 



Fig. 2. Block diagram of the dual-homodyne scheme for deducing the 
SQL for the resonant noise. 
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Fig. 3. Coherent State: Comparison of the smoothers for fix = 0.25 and 
= 0. 
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Fig. 4. Coherent State: Comparison of the smoothers for fii = 0.5 and 
= 0. 
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Fig. 5. Coherent State: Comparison of the smoothers for fi± = 0.75 and 
= 0. 



approaches —1, the performance of the robust smoother is 
superior to that of the Kalman smoother for all levels of fii. 
Also for fi\ — 0.25, the robust smoother is better than the 
SQL throughout the uncertainty window, while the Kalman 
smoother is not. 

Similarly, the estimation mean-square errors may be cal- 
culated, and plots created, for the robust smoother and the 



Kalman smoother as a function of the uncertain parameter 82 
with 81 — 0. Here, we choose different values for [12 and set 
jUi = 0. Figs. [6] and [7] show these plots for /12 = 0.2 and 0.5, 
respectively. In this case, the robust and Kalman smoothers 
behave almost identically for all levels of (12- 

VI. SQUEEZED STATE CASE 

We now consider a phase-squeezed beam as in Ref. [17], 
whose phase is modulated with the resonant noise process. 
The beam is then measured by homodyne detection using a 
local oscillator, the phase of which is adjusted according to 
the filtered estimate <pf(t). 

The normalized homodyne output current I(t) is given by 



where H 



I(t)dt~2\a\[(p(t)-<f>f(t)]dt + 



R s 



1 



-2r„ 



(28) 
(29) 



where |a| is the amplitude of the input phase-squeezed beam, 
and W(t) is a Wiener process arising from squeezed vacuum 
fluctuations. The parameter R sq is determined by the degree 
of squeezing (r m > 0) and anti-squeezing (r p > r m ) and by 



7 



The measurement model is, thus, [32]: 
9 = Hx + w, 



(30) 
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Fig. 6. Coherent State: Comparison of the smoothers for 112 = 0.2 and 
Ml = 0. 
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Eqs. (O to ( TT~8T > are then modified accordingly. 

VII. SQUEEZED STATE: COMPARISON FOR THE 
ROBUST AND KALMAN SMOOTHERS 

We use the Lyapunov method again, as described in sec- 
tion IV-AI to compute the estimation mean-square errors for 
the robust smoother and the Kalman smoother, as a function 
of 81 with 82 = 0. Here, we use the nominal values of the 
parameters, r rn = 0.36, r p = 0.59 and chosen values for [i\ 
with [12 = 0. These values were used to generate plots of the 
errors versus <5i to compare the performance of the robust 
smoother and the Kalman smoother for the uncertain system 
in comparison to the CSL. Due to the implicit dependence of 
R sq and <r 2 , we compute the smoothed mean-square error in 
each case by running several iterations until er 2 is obtained 
with an accuracy of 6 decimal places. The CSL value is 
obtained by designing, as described in section [ill] a different 
Kalman smoother for each value of the uncertain parameter 
in ( fT2b (note r m ,r p — and R sq — 1 here), and is given 
by ct 2 = P s (l, 1) from CHl- Figs. |8] |9] and [10] show the 
plots for [i\ — 0.3, 0.5 and 0.6, respectively. Note that as 8\ 
approaches —1, the performance of the robust smoother is 
better than that of the Kalman smoother for all levels of 

We can similarly plot the mean-square errors of the two 
smoothers compared to the CSL as a function of 82 with 
8\ = 0,/ii = for various values of /i2- Figs. ITTI and [T2l 
show these plots for /12 = 0.5 and 0.8, respectively. The 
Kalman and robust smoothers behave almost identically for 
all levels of fi2- 

VIII. CONCLUSION 

This work applies robust fixed-interval smoothing to ho- 
modyne phase estimation of coherent and squeezed states 
of light, when under the influence of a continuous-time 
resonant noise process. It presents a comparison of the robust 
smoother with the optimal smoother for the uncertain system 
in relation to the standard quantum limit for the coherent 
state case, and the coherent state limit for the squeezed state 
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Fig. 7. Coherent State: Comparison of the smoothers for fi2 = 0.5 and 
Ml = 0. 



Fig. 8. Squeezed State: Comparison of the smoothers for fii = 0.3 and 
M2 = 0. 
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Fig. 9. Squeezed State: Comparison of the smoothers for fii = 0.5 and 
= 0. 



Fig. 12. Squeezed State: Comparison of the smoothers for )i2 = 0.8 and 
Ml =0. 




Fig. 10. Squeezed State: Comparison of the smoothers for fii = 0.6 and 
,U2 = 0. 



case. These results may be demonstrated experimentally as 
further work. 
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